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Abstract. We present a model independent method of reconstructing the Lagrangian 
for the k-essence field driving the present acceleration of the universe. We consider 
the simplest k-essence model for which the potential is constant. Later we use three 
parametrizations for the Hubble parameter H(z), consistent with recent the SNla 
data, to yield the Lagrangian F. Our reconstruction program does not generate any 
physically realistic Lagrangian for models that allow phantom crossing, whereas models 
without phantom crossing, yield well behaved Lagrangians. 


PACS numbers: 98.80.-k, 98.80.Es, 98.70.Vc 

1. Introduction 

One of the most interesting observational discoveries in the past decade has been the 
evidence that the expansion of the universe is speeding up rather than slowing down. A 
currently accelerating universe is strongly favored by type la supernova data [I;. Recent 
observations of cosmic microwave background radiation (CMBR) by the Wilkinson 
Microwave Anisotropy Probe (WMAP) experiment p] and of the large scale structure 
by redshift surveys e.g the Sloan Digital Sky Survey (SDSS) j3] and the 2-Degree Field 
(2df) redshift survey |4J have further strengthened this result. This result together with 
the observational evidence that we are living in a low matter density (around 30% of the 
critical density) flat (D tota ; = 1) universe, suggests the existence of a relatively smooth 
component called “dark energy” which dominates the energy density of the universe 
and also has a large negative pressure. 

Although a simple cosmological constant A (a; = — 1) can serve the purpose of this 
dark energy, it faces a serious problem of fine tuning (for details, see recent review by 
Padanmabhan |5|). The alternative possibility called quintessence jUj, is a varying A 
model, where dark energy arises from a minimally coupled scalar field rolling over its 
potential. For a sufficiently flat potential, the potential energy term dominates and the 
scalar field can mimic A. 

Recently, an alternative possibility, that of an effective scalar field theory governed 
by a Lagrangian containing a non-canonical kinetic term (£ = —V{(j))F{ A), where 
A" = (1/2 )dfj,(j)d^(f>) has been proposed. Such a model can lead to a late time acceleration 
in the present universe and is named as “k-essence” [?]. It is worth noting that the 
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Generalized Chaplygin Gas model jSJ and the tachyon dark energy models |Jj may be 
seen as special cases of k-essence (e.g k-essence with a constant potential can mimic 
simple Chaplygin gas (CG) model [HI]). 

From the theoretical side, it is not only important to check whether dark energy 
is a constant or dynamical, but also if dynamical, it is equally important to know the 
Lagrangian for the underlying held. For a minimally coupled scalar held model, a 
number of studies have been done for reconstructing the potential for this scalar held 
in a model independent way by using observational data m- Recently Tsujikawa has 
studied the reconstruction of the potential function V ( 0 ) in the k-essence Lagrangian 
once one assumes a certain form for the kinetic energy function F(X) |12| . But until 
now, there is no attempt in the literature to reconstruct in a model independent way 
the function F(X) in the k-essence Lagrangian. 

In this work, we propose a algorithm (similar to that proposed by Simon et al. J32] 
for a minimally coupled scalar held model) to reconstruct F(X) in a model independent 
way. For present purposes, we restrict ourselves to the case where the potential V ( 0 ) is 
a constant. Using this algorithm and three specihc ansatzes for the Hubble parameter 
H(z) (consistent with the supernova observations), we reconstruct F(X) in the redshift 
range between z = 0 and z — 1.7 which is the current range for the supernova data. 


2. Cosmological reconstruction of F(X) 


In a spatially hat FRW universe, the luminosity distance Dl(z) and the coordinate 
distance r(z) to an object at redshift z are related as 

r= r*!m, (d 

Jt a(t) 1 + z 

where we have assumed the speed of light c = 1 and a 0 = 1, a{t) is the scale factor and 
quantities with subscript zero denote the values at present time (z — 0). This defines 
uniquely the Hubble parameter 


H(z) 


a 

d 

(D l (z)\] 

a 

dz 

U+* )\ 


( 2 ) 


This is a purely kinematic relation and does not depend on the choice of the 
matter part or the model of gravity. Hence knowing Dl{z) from observations, one 
can unambiguously determine the Hubble parameter H(z). 

Next, analogous to the inhationary Hubble-how parameter, we define the quantity 


_d log|e„| 

6n+1 dN ’ ^ 

for n > 0, where N = log[a(tj)/a(t)] is the number of e-foldings from some initial epoch 
ti and e 0 = jj-- One can now calculate C\ as 

dlog | e 0 1 E_ 

H 2 


ei = 


dN 


9 + 1, 


(4) 
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where q is the deceleration parameter. Similarly one can show 

d 


e 2 = -(1 + z) —[log(g + 1 )], 


( 5 ) 


which is related to the jerk parameter j = q+2q 2 +^l and also e 3 = — log(e 2 ) 

which is related to the snap parameter s [,Q| (These are similiar to the “Statehnder” 
parameters [T^) . We should emphasize that current supernova data can constrain up 
to the third derivative of the scale factor which means the deceleration parameter q and 
the jerk parameter j can be determined. 

Next, we define the energy density and pressure for a k-essence fluid. For a k-held 
cj) with a Lagrangian £ = —V((p)F(X), where X = > 0, the energy density 

and pressure are given by 

= V[F - 2 XF'} 

P*= - VF, 


( 6 ) 


where F' = jj^. One can now write Einstein’s equation for a flat FRW universe as 


3 H 2 — n 2 \p m + 


a 

a 


K r 
6 [pr 


+ P<t> + 3p<f ,\, 


( 7 ) 


where hi 2 = 87 tG and we neglect the radiation contribution as it is quite negligible around 
the present epoch and p m = Pbaryon + Pdm is the sum of the contributions coming from 
baryons and dark matter. 

In our subsequent calculations, we shall restrict ourselves to the simple k-essence 
models for which the potential V = Vo=constant. This case was first considered in PI 
and later discussed by a number of authors In this case one can 

show that 


= ° P m 7o\ 

1 2p m +[F-2XF’Y {) 

where we have assumed Vo = 1 without any loss of generality. Also, using Einstein’s 
equations and the expressions for p$ and p^, one can write 


-2 XF 


tt2 

F — —y [3 — 2 ei] 


K 


Using eqns ( 8 ) and (9), one can now write 

C\H~ p m 

' ~ 2 ' 


-XF = 


K 


(9) 


( 10 ) 


In our subsequent calculations, we shall use eqns (9) and (10) for reconstructing 
F(X). Next using the relation e\ = one can write eqn (9) as 


F = fr[ 3 -% + 1 )]’ 

which at the present epoch becomes 
H 2 

To = —y[3 — 2(go + 1)]- 


( 11 ) 


( 12 ) 
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Figure 1. Reconstructed F for Ansatz 1 and 2. The solid line is for ansatz 1 and the 
dashed line is for ansatz 2. 


Hence knowing H(z) from observation, one can reconstruct F(z). Also using eqn 
(9) one can determine F' (z)\q (here prime denotes derivative w.r.t the argument z): 

F ( z ) = k 2 l + z ^ _ 261 + ^ ’ ( 13 ) 

which gives F(z) at present: 

F (z)\o = -p^[3(<?o + 1) ~ 2 (g 0 + l) 2 - (^r)lo]- (14) 

Similarly one can determine F"(z),F"\z) and so on at present epoch. One can always 
make a Taylor expansion of F(z) around the present epoch z — 0: 

F(z) = Fq + (F |o )z + i^(F |q )z 2 + .... (15) 

Hence it is always possible to determine the approximate behavior of F(z) around 
the present epoch by knowing H(z), q(z), q'(z), q"{z) etc at present (z = 0) from 
supernova observations even if one does not know the full evolution of H{z). Current 
supernova data can only determine up to q (z = 0) [21, hence in this way one can 
determine F{z) around the present epoch up to only linear order. But as future 
observations start probing the universe at higher redshifts, one can determine the higher 
derivatives of q(z), and hence the higher order behavior of F(z). This provides a nice 
and useful method of determining the behavior of F(z) around the present epoch from 
future distance-redshift observations. 
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But our goal is to reconstruct F(X), which is not what eqn (11) or (15) provides. 
For this, we consider the energy conservation equation for the K-Field: 

(F'(X) + 2 XF"(X))4> + 3 HF\X)<j) = 0 (16) 

which gives a first integral for F(X) as [Rh IT71 21] 

XF" 2 (X) = c( 1 + z) 6 , (17) 


where c is an arbitrary integration constant. One can now use eqn (10) to write c as 


c = 


H o 


«r 




rO m o ^i( z 0) 


(18) 


where X 0 is the present day value of X and h2 m ,o is the present day density parameter 
for matter. Using this expression for c, one can now write 

ei (l + ^)- 31 


X 

TO 


= E\z ) 


3 ^mO 
2 E 2 (z) 


[f^mO ~ e l( z ~ 0)] 


(19) 


where E(z) = jj-. Eqns (11) and (19) constitute the main equations for our 
reconstruction procedure. Eqn. (11) gives F as a function of z and Eqn. (19) gives 
X as a function of Z. Hence together, they determine F(X). This is important as we 
want to know F(X) not F(z). 

As mentioned above, our method of successful reconstruction of F(X) depends on a 
sufficiently versatile fitting function for H(z). For our purpose we shall use the following 
flexible and model independent ansatzs: 


Ansatz 1: 


H(z ) — H o[n m0 (l + ^) 3 + A) + A\ (1 + z) 

+ A 2 {l + z) 2 ] 1 ' 2 . (20) 

This has been proposed by Alam et al. [23, [24J. For a flat universe, the model has three 
independent parameters as Vl m0 +A 0 +A 1 + A 2 = 1. For A 1 = A 2 = 0 the ansatz becomes 
an exact ACDM model. For A 0 = A 2 = 0 and for A 0 = Ai = 0, it gives a dark energy 
model with u = —2/3 and to = —1/3 respectively. It has also been shown by Alam et 
al. |22] that even for models where this ansatz is not exact, like tracker quintessence 
models, supergravity models or k-essence models like the Chaplygin gas, one can recover 
the luminosity distance to within 0.5% accuracy using this ansatz in the redshift range 
that is relevant for supernova observations. Hence one can trust this ansatz to recover 
cosmological quantities to a high degree of accuracy for a large class of models. Best fit 
values for the parameters Ai and A 2 for fl m = 0.3 are Ai = —4.16 and A 2 = 1.67, for 
the SNla Gold dataset [23] and will be used for subsequent reconstruction. 
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Figure 2. Variation of X/Xq (solid line), F (dashed line) and equation of state for 
k-essence field (dot-dashed line) for Ansatz 1 


Ansatz 2: 

w{z) = w 0 + ( 21 ) 

1 + z 

where Wq and W\ are constants. This ansatz was first discussed by Chevallier and 
Polarski [2BJ and later studied more elaborately by Linder [25|. For low redshifts, 
the ansatz reduces to the conventional first order expansion to the equation of state 
w = Wo + W\z. It is also well behaved and bounded for high redshifts and can mimic 
several scalar field equations of state. For example, for the supergravity model, the 
ansatz can mimic the distance-redshift behavior to 0.2% up to redshift 0 ~ 1100. The 
best fit values for this model while fitting with the SNla gold dataset are wq = —1.58 
and w\ = 3.29 for Q m = 0.3 233- We shall use these values for the parameters for the 
reconstruction. 

Ansatz 3: 

H(z) = H 0 [fl m0 (l + z) 3 

+(1 - fl m0 )(A s + (1 - A a )(l + ^)3(i+«))i/(i+«)]i/2 (22) 

where A s and a are constant. This ansatz for the Hubble parameter H{z) can be 
obtained by modeling the dark energy as Generalized Chaplygin gas with an equation 
of state p c h = —Jr [§j. By choosing different ranges for the parameters A s and a, 
this ansatz can behave as standard dark energy model with asymptotic de-Sitter phase, 
early phantom model, late phantom phantom model as well as transient model where 
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Figure 3. Variation of X/Xq (solid line), F (dashed line) and equation of state for 
k-essence field (dot-dashed line) for Ansatz 2 


the present acceleration of the universe is only temporary as it again enters the dust- 
dominated decelerating phase in future. Recently it has been shown that current SNla 
gold data can allow all four behaviors of this model with suitable choice of the matter 
density parameter Q m0 [27,. For our analysis, we have considered four combinations of 
(A s ,a): (0.95, 0.2), (0.85, -1.2), (1.1 -0.8) and (1.5, -1.7), which result in the standard, 
transient, early phantom and late phantom phases respectively for = 0.2, 0.25, 0.3 
and 0.35. All these combinations for (A,, a) are within the la confidence contour for 
SNla Gold dataset 27. 

Results: 

In Fig 1, we have plotted the reconstructed F as a function of X/X 0 for ansatz 1 and 
2. For both of these ansatzes, the reconstructed F(X) is a double-valued function of 
X. This behaviour of F(X) is not what one should expect for any physically realistic 
Lagrangian. To see the reason behind this pathological behaviour for F(X), one has 
to carefully study the equations. First of all, one can see both of these ansatzes allow 
the crossing of the phantom divide (a; = —1) [2B3- Also X(z) has a minimum in both 
cases with X m i n = 0 and the minimum occurs precisely when the equation of state 
crosses the phantom divide uj — — 1 (see Fig 2 and 3). On the other hand, F is a 
monotonically decreasing function with z within the relevant redshift range. This forces 
F to be a double-valued function of A" when F doubles back on itself at X = 0. Hence 
one might suspect that this unphysical behaviour of F(X) is due to the crossing of 
phantom divide. To see that this is actually the reason, let us consider the equation of 
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Figure 4. Reconstructed F(X) for anstaz 3 for four different cases (See text) 


state of the k-essence fluid: 


u = 


nx) 

2 XF'(X) - F(xy 


(23) 


In order to have u = —1, one should have X = 0, or F'(X) = 0, or both. But looking at 
the equation of motion (eqn 17), one concludes that it is not possible (except at future 
infinity z = —1), unless X = 0 and F'{X) —> oo as A=, so that the l.h.s of eqn. 17 
remains finite. Also, to cross the boundary u = —1, F' {X) should change its sign at 
X = 0. These two conditions (e.g F\X) diverges and changes sign at X = 0) together 
single out the sole possibility of F(X) being a double-valued function and crossing 
back on itself at A" = 0, the point of phantom crossing. This rules out any physically 
realistic Lagrangian F(X) for a k-essence held which allows phantom crossing. This is 
consistent with earlier results by Vikman [22] and Caldwell and Doran EH who have 
made a more elaborate study for phantom crossing with a k-essence held and arrived at 
similar conclusions, (see also the more recent investigation fcTTl.) 

But one can check that restricting the equation of state from crossing the phantom 
divide can result in a well-behaved reconstructed F(X). This happens for our ansatz 
3 which although perfectly consistent with the recent SNla data, never crosses the 
phantom divide. In this case the reconstructed F(X) is shown in Fig 4 for four different 
models and it can be seen that F(X) is a perfectly well-behaved single valued function 
for all four cases. It is also interesting to note that for standard and transient models 
(non-phantom cases), F is a monotonically decresing function of X whereas for early 
and late phantom models, it is an increasing function of A" (see Fig. 4). This is because 
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Figure 5. Variation of X(z) (solid line) and F(z) (dashed line) for four different cases 
(See Text) for ansatz 3 


for X > 0, the sign of is related to the value of to (see Eqn.(23)) and for phantom 
models, F'(X ) > 0, whereas for non-phantom case, F'(X ) < 0 as F(X) is always 
positive for all four cases in the relevant redshift range (see Fig 5). 

The behaviors of X(z) and F(z) for this ansatz are shown in Fig 5. 

3. Conclusion 

In this work, we have proposed a method for reconstructing the k-essence Lagrangian 
for constant potential. In doing so, we have defined the parameters e n , which are similar 
to the horizon-flow parameters for inflation. Using these parameters, we have related 
F(z) and F'(z) to the directly observable quantities like Q m , q(z), H(z). Equations (11) 
and (19) are the core of our reconstruction program. By knowing H(z) from distance- 
redshift measurements like SNIa observations, one can reconstruct F as a function of 
redshift z and X, the kinetic energy. Later, we have used three well known ansatz for 
H(z ), in order to reconstruct F for z = 0 to z = 1.7 which is the current redshift range 
for SNIa data. We should emphasize that our reconstruction program only involves that 
portion of F(X) over which the field evolves to give the required H(z) and is insensitive 
to certain ranges in X. 

Our study shows that in order to have a physically realistic k-essence Lagrangian 
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F(X), the crossing of the phantom boundary u = — 1 is forbidden. For ansatzes 1 and 
2, which allow the phantom crossing, the reconstructed F ) s are physically unrealistic 
double-valued functions of X. This result is consistent with earlier investigations with 
similar conclusion HU imi EH- For ansatz 3, which is also consistent with SNla data 
but does not allow phantom crossing, the reconstructed F(X) is a well-behaved, single 
valued function. The result also shows that one has to be careful while choosing different 
ansatzes for H(z) or u(z) for dark energy to fit different observational data. While some 
of them are perfectly consistent with observation, it may not be possible to relate them 
to any physically realistic theory. 
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